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Based on a data sample of 106 M ip' events collected with the BESIII detector, the decays 
ifr' — > 7XcO,2,XcO,2 — > 77 are studied to determine the two-photon widths of the Xc0,2 states. The 
two-photon decay branching fractions are determined to be B(xco — > 77) = (2.24 ± 0.19 ± 0.12 ± 
0.08) x 10~ 4 and B(xc2 -+77) = (3.21±0.18±0.17±0.13) x 10~ 4 . From these, the two-photon widths 
are determined to be r 77 (xeo) = (2.33±0.20±0.13±0.17) keV, T 77 (x C 2) = (0.63±0.04±0.04±0.04) 
keV, and TZ = r 77 (xc2)/r 77 (xco) = 0.271 ± 0.029 ± 0.013 ± 0.027, where the uncertainties are 
statistical, systematic, and those from the PDG B(ip' — > 7Xc0,2) and T(xc0,2) errors, respectively. 
The ratio of the two-photon widths for helicity A — and helicity A = 2 components in the decay 
Xc2 -> 77 is measured for the first time to be f 0/2 = r^ (xc2)/r^= 2 (x C 2) = 0.00 ± 0.02 ± 0.02. 

PACS numbers: 13.20. Gd, 13.25.Gv, 14.40.Pq, 12.38.Qk 



I. INTRODUCTION 

Charmonium physics is in the boundary domain be- 
tween perturbative and nonperturbative quantum chro- 
modynamics (QCD). Notably the two-photon decays of 
P-wave charmonia are helpful for better understanding 
the nature of interquark forces and decay mechanisms [l[ . 
In particular, the decays of Xco.2 — ► 11 offer the closest 
parallel between quantum electrodynamics (QED) and 
QCD, being completely analogous to the decays of the 
corresponding triplet states of positronium. In the low- 
est order, for both positronium and charmonium the ra- 
tio of the two-photon decays TZ^ = ^(^^77) = 4/15 w 
0.27 Any discrepancy from this simple lowest or- 
der prediction can arise due to QCD radiative correc- 



tions and relativistic corrections, and the measurement 
of TZ provides useful information on these effects. The 
decay of \ci — > 11 is forbidden by the Landau- Yang the- 
orem 0. Theoretical predictions on the decay rates are 
obtained using a non- relativistic approximation po- 
tential model @ , relativistic quark model @, H| , nonrel- 
ativistic QCD factorization framework effective La- 
grangian [l(|, as well as lattice calculations [ll|. The 
predictions for the ratio TZ = r^^fco) cover a wide range 
values between 0.09 and 0.36 [H, Iff- Precision measure- 
ments of these quantities will guide the development of 
QCD theory. 

The two-photon decay widths of \cj have been mea- 
sured by many experiments [l2j |. Using the reactions 
ip' — > iXcj, the CLEO-c experiment reported results for 
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r 77 (y £ j) measured in the decay of XcJ into two pho- 
tons 0: 

r 77 ( Xc0 ) = (2.36 ± 0.35 ± 0.22) keV, (1) 
r 77 ( Xc2 ) = (0.66 ± 0.07 ± 0.06) keV, 

with uncertainties that are dominated by the statistical 
errors. BESIII has collected 106 million ip' events, a data 
sample that is about four times of that of CLEO-c, allow- 
ing for more precise measurements of these quantities. 

There are two independent helicity amplitudes, the 
helicity-two amplitude (A = 2) and the helicity-zero 
(A = 0) amplitude, that contribute to Xc2 — > 77 de- 
cay 0, where A is the difference in the helicity val- 
ues of the two photons. The ratio of the two-photon 
partial widths for the two helicity components, / /2 — 
r 77 (x C 2)/r^ 2 (Y C 2) in the decay \ci -> 77: is predicted 
to be about 0.5% @; a measurement of this ratio can be 
used to test the QCD prediction. 

In this paper, (1.06±0.04) x 10 s ip' events accumulated 
in BESIII are used to study the process ip' — > 7iXeO,2, 
XcQ.2 — > 7273 and measure the two-photon decay widths, 
r 77 (x c o) and r 77 (x C 2)- We also determine the ratio 1Z, 
where many of the systematic uncertainties cancel in the 
ratio of the two simultaneous measurements. The ra- 
tio of the helicity-zero component relative to helicity-two 
component, /o/2> is also reported for the first time. 

II. THE BESIII EXPERIMENT AND DATA SET 

This analysis is based on a 156.4 pb" 1 of ip' data cor- 
responding to (1.06 ±0.04) x 10 8 ip' events pj| collected 
with the BESIII detector [lj| operating at the BEPCII 
Collider [l6j]. In addition, an off- resonance sample of 44.1 
pb _1 taken at y 7 ! = 3.65 GeV is used for the study of 
continuum backgrounds. 

BEPCII/BESIII is a major upgrade of the BE- 
SII experiment at the BEPC accelerator [T3] for studies 
of hadron spectroscopy and r-charm physics 18]. The 
design peak luminosity of the double-ring e + e~ collider, 
BEPCII, is 10 33 cm" 2 s" 1 at a beam current of 0.93 A. 
The BESIII detector has a geometrical acceptance of 93% 
of 47r and consists of four main components: (1) a small- 
celled, helium-based main draft chamber (MDC) with 43 
layers. The average single wire resolution is 135 fim, and 
the momentum resolution for 1 GeV/c charged particles 
in a 1 T magnetic field is 0.5%; (2) an electromagnetic 
calorimeter (EMC) made of 6240 Csl (Tl) crystals ar- 
ranged in a cylindrical shape (barrel) plus two end-caps. 
For 1.0 GeV photons, the energy resolution is 2.5% in 
the barrel and 5% in the end-caps, and the position res- 
olution is 6 mm in the barrel and 9 mm in the end-caps; 
(3) a time-of-flight system (TOF) for particle identifica- 
tion composed of a barrel part made of two layers with 
88 pieces of 5 cm thick, 2.4 m long plastic scintillators in 
each layer, and two end-caps with 96 fan-shaped, 5 cm 
thick, plastic scintillators in each end-cap. The time res- 
olution is 80 ps in the barrel, and 110 ps in the end-caps, 



corresponding to a 2a K/7r separation for momenta up to 
about 1.0 GeV/c; (4) a muon chamber system made of 
1000 m 2 of resistive plate chambers arranged in 9 layers 
in the barrel and 8 layers in the end-caps and incorpo- 
rated in the return iron of the super-conducting magnet. 
The position resolution is about 2 cm. 

The optimization of the event selection and the estima- 
tion of physics backgrounds are performed using Monte 
Carlo (MC) simulated data samples. The GEANT4-based 
simulation software BOOST [19| includes the geometric 
and material description of the BESIII detectors, the de- 
tector response and digitization models, as well as the 
tracking of the detector running conditions and perfor- 
mance. The production of the ip' resonance is simu- 
lated by the Monte Carlo event generator KKMC 1201 : the 
known decay modes are generated by evtgen [21| with 
branching ratios set at PDG [l2[ world average values, 
and by lundcharm [22j for the remaining unknown de- 
cays. The analysis is performed in the framework of the 
BESIII offline software system [23[ which takes care of 
the detector calibration, event reconstruction and data 
storage. 



III. DATA ANALYSIS 

Electromagnetic showers are reconstructed from clus- 
ters of energy deposits in the EMC crystals. The en- 
ergy deposited in nearby TOF counters is included to 
improve the reconstruction efficiency and energy reso- 
lution. Showers identified as photon candidates are re- 
quired to satisfy fiducial and shower-quality criteria. A 
photon candidate is a shower detected in the EMC with 
a total energy deposit greater than 25 MeV and with an 
angle 9 with respect to the e + beam direction in the range 
| cos#| < 0.75. This requirement is used to suppress con- 
tinuum background e + e~ — > 77(7), where the two ener- 
getic photons are mostly distributed in the forward and 
backward regions. We restrict the analysis to events that 
have no detected charged particles. The average event 
vertex of each run is assumed as the origin for the se- 
lected candidates. For ip' — > 7iXc0.2, Xc0.2 — > 7273 analy- 
sis, events are required to have three photon candidates, 
among which the smallest energy photon is selected as 
the radiated photon 71 and the second-largest and the 
largest energy photons are defined as 7273 from Xc0,2 de- 
cays. An energy-momentum conservation constraint 4C- 
fit is performed, and events with x 2 < 80 are retained in 
the final selection. The energy spectrum of the radiated 
photons is shown in Fig. [TJ where enhancements due to 
the XcO and Xc2 over substantial backgrounds are clearly 
observed. 

To determine signal efficiencies 100K signal MC event 
samples are generated for the XcO and the Xc2, with PDG 
values for the masses and widths [l2j ■ The radiative tran- 
sition ip(2S) — > 7iXc0 is generated using a (1 -I- cos 2 #) 
distribution, where 9 is the radiative photon angle rela- 
tive to the positron beam direction, in accordance with 
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TABLE I: Expected numbers of background events peaking 
at the Xcj signal regions from MC simulations. The errors 
are the uncertainties from these measured branching frac- 
tions [3. 



Decay modes 




n XcO 




V>' ->■ 7XcO,XcO - 


+ 7rV 


25.4 ±2.2 


0.0 ±0.0 


i>' 7XcO,XcO - 


-»■ 7777 


0.4 ±0.1 


0.0 ±0.0 


i>' ->■ 7Xc2,Xc2 - 


■> 7T 7T 


0.0 ±0.0 


7.7 ±0.7 


i>' -> 7Xc2,Xc2 - 




0.0 ±0.0 


0.1 ±0.1 


Sum 




25.8 ±2.2 


7.8 ±0.7 



expectations for pure El transitions. The X c0 ~~ ^ 7273 
decays are generated using a uniform angular distribu- 
tion. Although the radiative transition ip(2S) — > 7iXc2 
is dominantly pure El [24|,[25j], there is some recent exper- 
imental evidence that the decay has contributions from 
higher-order multipoles [26J . The full angular amplitudes 
for ip' — > 7ix C 2 are discussed in association with Eq. ([5]) 
in Section|Vl Furthermore, the 7273 photons in the decay 
Xc2 7273 are expected to be mostly in a pure helicity- 
two state; the ratio of the partial two-photon widths for 
the helicity-zero and helicity-two amplitudes is predicted 
to be less than 0.5% @- Thus the signal MC for the de- 
cay ip' -> 7iXc2, Xc2 -> 7273 is generated with 7273 in a 
helicity-two state as described in Section [V] 

The energy resolutions determined by the MC sim- 
ulations are <r(E 7l ) = 6.74 ± 0.29 MeV for Xco and 
a(E 71 ) = 3.91 ± 0.09 MeV for X c2- The efficiencies de- 
termined from MC simulations for the Xco and \ci are 
e(x c o) = (35.4±0.06)%ande(x c2 ) = (38.0±0.07)%. The 
difference between e(xco) and e(Xc2) is due primarily to 
the different angular distributions. 

The dominant non-peaking background that is ap- 
parent in the spectrum in Fig. [T] is from continuum 
e + e~ — > 77(7) processes. It is determined from MC 
simulations that contributions to the background due to 
radiative decays to the 77, rj , and 37 decays of ip' are non- 
peaking, spread over the full range of E- fl , and negligible. 
Therefore, they do not change the shape of the dominant 
continuum background. In addition we use MC simu- 
lations to investigate possible sources of peaking back- 
grounds. These are found to come from Xc0.c2 7r°7r° 
and rjij decays and Tr°(r)) — > 77, where two of the 7s 
have low momentum and are not detected or are outside 
of the fiducial volume of this analysis. We generate at 
least 100K events of each type to determine the efficien- 
cies for the peaking backgrounds, and use the efficiencies 
and branching fractions measured by BESIII [l4| to de- 
termine the numbers of peaking background events listed 
in Table [J 



IV. MEASUREMENT OF BRANCHING 
FRACTIONS AND TWO-PHOTON WIDTHS 

An unbinncd maximum likelihood (ML) fit is done 
to the E yi spectrum as shown in Fig. [1] The shape 
of the large nonpeaking background in the spectrum is 
determined with the 44.1 pb^ 1 of off-ip' data taken at 
y/s = 3.65 GeV, as well as the 921.8 pb" 1 of 0(3770) 
data taken at \/s — 3.773 GeV. As is evident in Fig. [5J 
the o&-ip' data are in good agreement with the high statis- 
tics V>(3770) data, for which transitions to either the x r n 
or X c2 states are expected to be less than 8 events (12] ]. 
We also generate e + e~ —> 77(7) MC events using the 
Babayaga QED event generator (2tJ and confirm that 
the shapes from the 3.65 GeV and 3.773 GeV samples 
are consistent with being due to the QED process. The 
E 7l distribution for the V(3770) data is fitted with the 
data-driven function: 

f bg (E 71 ) = po + pi x E 71 ± p 2 x (E 71 ) a , (2) 

where po, p\, pi and a are parameters which are obtained 
in a fit to the '0(3770) data in Fig. [2] In the nominal fit to 
the data, the background shape is fixed to Eq. @ , but 
its normalization is allowed to float. The shapes of the 
Xco and Xc2 resonances used in the fit are extracted from 
a nearly background-free tp' — !• 7iXc0.2, Xc0,2 — > K + K~ 
sample shown in Fig. [3] The purity of the sample is larger 
than 99.2%. The shapes of the signal peaks in the E*yi 
spectrum arc fixed to the smoothed-histograms of the 
tp' — > 7iXcO,2, Xc0,2 — > K + K~ sample, and the yields are 
allowed to float. The estimated numbers of peaking back- 
ground events from Xc o,c2 7r°7r° and rjr] that contribute 
to the Xco and Xc2 signals are 25.8 and 7.8 events, respec- 
tively, as listed in Table |U They are subtracted from the 
fitted yields, and after this subtraction, the signal yields 
are N( Xc o) = 813 ± 63 and N{ X c2) = 1131 ± 66. The 
product branching fractions are determined from the re- 
lation 

B(V' 7V./) x B( Xc j 77) = 7^%4 , (3) 

e (XcJ) x iV^/ 

where N^i is the total number of in the data sample. 
The measured product branching fractions are listed in 
Table ILT1 We use the PDG average values, 

£(</>'-> 7Xco) - (9.68 ± 0.31) x 10- 2 , (4) 

rfeo) - (10.4±0.6)MeV, 

W^7Xc 2 ) = (8.75 ± 0.35) x 10- 2 , 

r(x c2 ) = (1.97±0.11)MeV, 

to determine B(xc0,2 — ► 77), r 77 (x c o,2) and 1Z. These 
are also listed in Table |TTJ 

Several sources of systematic uncertainties in the mea- 
surement of the branching fractions are considered, in- 
cluding: uncertainties on the photon detection and re- 
construction; the number of tp' decays in the data sample; 
the kinematic fitting; the fitting procedure and peaking 
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FIG. 3: The E 71 spectrum for the radiative photon in the 
samples ip' 7iXcO,2, Xc0,2 -t K + K~ . 



FIG. 1: Upper plot: the fitted i? 71 spectrum for the rp' data 
sample. The expected positions of E~ tl from XcO, Xci, Xc2 are 
indicated by arrows. The dashed curve shows the background 
line shape fixed to the shape in Fig. [5] Lower plot: the num- 
ber of standard deviations, x, of data points from the fitted 
curves. 




0.10 0.15 0.20 0.25 0.30 0.35 0.40 

E ( Yl ) (GeV) 

FIG. 2: The background E 11 spectrum. The points are from 
the off-ip' data. The curve is from a fit to the ^(3770) data. 



background subtraction. Table |TTT] lists a summary of 
all sources of systematic uncertainties. Most systematic 
uncertainties are determined from comparisons of spe- 
cial clean, high statistics samples with results from MC 
simulations. 

The number of ip' events, N^> , used in this analysis 
is determined from the number of inclusive hadronic ip' 



TABLE II: Results of the present measurements. The first 
error is statistical, second is systematic, and third is due to 
the PDG values used. The common systematic errors have 
been removed in determining 1Z. B\ = B(tp' — > 7X00,2), B2 = 
B(xco,2 ->■ 77). L 77 = r 77 (x c o,2 ->■ 77). 



Quantity 



XcO 



Xc2 



BixB 2 x 10 s 
B 2 x I0 4 



(keV) 



n 



2.17±0.17±0.12 2.81±0.17±0.15 
2.24±0.I9±0.12±0.08 3.21±0.18±0.17±0.13 
2.33±0.20±0.13±0.I7 0.63±0.04±0.04±0.04 
0.271±0.029±0.013±0.027 



TABLE III: Summary of systematical uncertainties of the 
branching fraction measurements. Asterisks denote the sys- 
tematic sources common to both XcO and Xc2- 



Source of Systematic Uncertainty 


XcO 


Xc2 


Number of tp'* 


4.0% 


4.0% 


Neutral trigger efficiency* 


0.1% 


0.1% 


Photon detection * 


1.5% 


1.5% 


Kinematic fit * 


1.0% 


1.0% 


Resonance fitting 


3.2% 


2.9% 


Peaking background 


0.3% 


0.1% 


Helicity 2 assumption 




0.4% 


Sum in quadrature 


5.5% 


5.3% 



decays following the procedure described in detail in [14) . 
The result is N^> = (1.06 ± 0.04) x 10 8 , where the error 
is systematic. 

Three photons in the final states include a soft pho- 
ton 71 from the radiative transition and two energetic 
photons 7273 from Xc0,2 decays. The photon detection 
efficiency and its uncertainty for low energy photons 
are studied using three different methods described in 
Ref. [28| . On average, the efficiency difference between 
data and MC simulation is less than 1% [28]. The mo- 
menta of the two energetic photons are more than 1.5 
GeV/c. The systematic uncertainty due to the recon- 
struction of two energetic photons is determined to be 
0.25% per photon as described in Ref. 29]. The total un- 
certainty associated with the reconstruction of the three 
photons is 1.5%. 

The uncertainty due to the kinematic fit is estimated 
using a sample of e + e~ — > 77(7), which has the same 
event topology as the signal. We select the sample by 
using off--0' data taken at \fs = 3.65 GeV to determine 
the efficiency difference between data and MC for the 
requirement of \\c < ^0 in the 4C-fit. The uncertainty 
due to kinematic fitting determined in this way is 1%. 

Since the signal shapes are obtained from ip' — > 7X c o,2, 
XcO, 2 — > K + K~ events in the data, the uncertainty due 
to the signal shape is negligible. The shape of the contin- 
uum background is parameterized using the data-driven 
function in Eq. @; the parameters obtained in the fit- 
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ting to oW-tp' data sample are fixed in the nominal fitting 
to data. The systematic uncertainty due to the choice 
of parametrization for the background shape is estimated 
by varying the fitting range and the order of polynomial 
in our data-driven function. We find relative changes on 
the XcO and %c2 signal yields of 3.2% and 2.9%, respec- 
tively, which are taken as the uncertainties due to the 
resonance fitting. 

The expected numbers of peaking background events 
from Xci.2 — > 7T 7r° and Xc0.2 - > VV decays summa- 
rized in Table U use BESIII measurements for S(Xci,2 — ^ 
ir° it / rjr]) [l4j]. The uncertainties on the 7r°7r / 7777 contri- 
butions are estimated to be 0.3% and 0.1% for XcO and 
Xc2, respectively. The systematic uncertainties due to 
the trigger efficiency in these neutral channels are esti- 
mated to be < 0.1%, based on cross-checks using differ- 
ent trigger conditions [H], Ho|. We have assumed pure 
helicity-two decay of \ci ~ ► 77- In a relativistic calcu- 
lation, Barnes [5| predicted the helicity-zero component 
to be 0.5%. In section [Vj the ratio of the two photon 
widths for the helicity-zero and helicity-two amplitudes 
is measured to be 0.00 ± 0.02 ± 0.02. To be conserva- 
tive, we determine the change in our Xci result when a 



helicity-zero component of 3% is included, corresponding 
to an upper limit at 90% confidence level from the mea- 
surement in this paper, to be 0.4%, and use that as the 
helicity-state-associated systematic error. 

All sources of systematic errors are listed in Table IIIII 
We assume that all systematical uncertainties are inde- 
pendent and add them in quadrature to obtain the to- 
tal systematical uncertainty. For the measurements of 
B{XcO,2 — > 77), the uncertainty due to the ip' —> 7Xc0, 2 
branching fractions is kept separate and quoted as a sec- 
ond systematic uncertainty. 



V. HELICITY AMPLITUDE ANALYSIS FOR 

Xc2 ->■ 77 

In Xc2 — > 77 decay, the final state is a superposition 
of helicity-zero (A = 0) and helicity-two (A = 2) com- 
ponents, where A is the difference in the helicity values 
of the two photons. The formulae for the helicity am- 
plitudes in ip' — > 7Xc2,Xc2 — > 77, which include higher- 
order multipole amplitudes, are: 



3 •}.,, , 2„n.4„ , „ H . H „ ■ 2 „ n 3\/2 _ 



y (1 + cos 0i)sin 82 + 3x sin 0i sin 282 — :r j/ sin 20i sin 82 sin 202 cos <p2 



W2(0 1 ,8 2 ,<h) = fo/2 

+ \/3x sin 28j sin 28 2 (3 cos 2 8 2 - 1) cos <f> 2 + -My sin 2 81 sin 2 2 (3 cos 2 8 2 - 1) cos 202 + (1 + cos 2 6»i ) (3 cos 2 2 - l) 2 



1 \f2 

-2/ 2 (l +cos 2 6»i)(l + 6 cos 2 82 + cos 4 82) + 2x 2 sin 2 0i(l + cos 2 2 )sin 2 8 2 + ^-xysin28 1 sin20 2 (3 + cos 2 2 ) cos<fe 



x sin 20i sin 2 02 sin 20 2 cos <j>2 + — — y sin 2 8\ (1 — cos 4 2 ) cos 2^> 2 H — (1 + cos 2 0i) sin 4 02 



(5) 



q 

where x = Ai/A , y = A 2 /A , and A ,i, 2 are the x c2 fl2 = -(1 + cos 2 #1) sin 4 6» 2 , (7) 

helicity 0, 1, 2 amplitudes, respectively, d\ is the polar 2 
angle of the radiative photon, where the electron beam 

is defined as the z axis in the e + e~ center-of-mass frame, r- 

and 6» 2 and 4> 2 are the polar angle and azimuthal angle of fl3 = 1— sm 20 1 sin 2 6 2 sin 26 2 cos (j) 2 , (8) 

one of the photons from Xc2 decay in the Xc2 rest frame, 2 
relative to the radiative photon direction as polar axis; 

62 = is defined by the electron beam direction. Thefac- n: . nn . „„ ,„ 2n 1N , lc . s 

tor f Q/2 = |F | 2 /|f 2 | 2 = r^(x c2 )/r^ 2 (x c2 ) is the ra- &4 = VSsn^ sm20 2 (3cos 9 2 - 1 )cosfe (9) 

tio of partial two-photon widths for the helicity-zero and 

helicity-two components, where Fn (F 2 ) is the helicity- fc ■ 2 n ■ 2 n m 2 n ,1 nj. /m\ 

fl s 1 • a , t? +1, a 5 = v6srn 6»i sin 6> 2 (3cos z 9 2 - 1) cos2</> 2 , (10) 

zero (two) amplitude m the decay Xc2 ~^ 77- rurtner 

information on the formulae for the helicity amplitudes 

can be found in Ref. [p. a 6 = (1 + cos 2 X )(3 cos 2 9 2 - l) 2 , (11) 

An unbinned ML fit to the angular distribution is per- 
formed to determine x, y and f 0/2 values. We define o^fl h j_ ™c 2 a \ a \^ft <-\o\ 
twelve factors j2fj: ° 7 = 2 Sm + C ° b 02 > Sm ^ (12) 

a x = 3sin 2 1 sin 2 20 2 , (6) a 8 = J(l + cos 2 ^(l + 6 cos 2 2 + cos 4 2 ), (13) 
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v2 .„„.„„/„, 2/i \ i / j \ used for the simulation of all the angular distributions: 

— sm 20i sin 20 2 (3 + cos 02 J cos 02, (14) 



aio = — — sin 20i sin 2 02 sin 202 cos 02, (15) 



y~L i a« (*) 

~ an= N » w = 1 '--- > 12 



(18) 



V^6 

an = -— sin 2 0i (1 — cos 4 02) cos 202, 



ai2 = -(1 + cos 2 0i) sin 4 2 . 



(16) 

where A^ is the number of events after all selections from 
phase space MC samples. Since a n is calculated with 
phase space MC events after selection, it naturally ac- 
(17) counts for the detector acceptance effects. 



The mean values of ai, • • • , ai2 can be determined with 



The normalized probability-density function is written 



V*' ~~ * 7Xc2, Xc2 — > 77 MC events, where phase space is as: 
I 



f(x,y,fo/2) = 



W 2 (di,6 2 ,<t> 2 \x,y,f 0/2 ) 



fo/ 2 {aix 2 + a 2 y 2 + a 3 xy + a^x + a 5 y + a 6 ) + (a 7 x 2 + a$y 2 + a 9 xy + a w x + any + ai 2 ) 



(19) 



r 



A total log-likelihood function is constructed as: ln£ — 
Y^i=i m fi( x i 2/, fo/ 2 )y where the sum is over all the events 
in the signal region (here the signal region is defined as 
0.09 < E 71 < 0.15 GeV). The log-likelihood function for 
the signal is given by In C s — In C — In Ci, , in which In is 
the normalized sum of logarithmic likelihood values from 
background events and is calculated using the events in 
the sidebands, which are defined in the ranges (0.07, 0.08) 
GeV (lower sideband) and (0.16,0.20) GeV (higher side- 
band) in the the E lx spectrum. By maximizing the loga- 
rithm of the likelihood function ln£ s , the best values of 
x, y and /0/2 are determined. Before fitting to the data, 
input and output checks were done using MC samples, 
and the results used to validate the fitting procedure. 

BESIII has determined x and y to be x — 1.55 ± 
0.05(stat.) ± 0.07(syst.) and y = 2.10 ± (stat.)0.07 ± 
0.05(syst.) [Ifl using the decays tp' — > "fXc2, Xc2 — > 
tt + n~ / K + K~ . Therefore, in the nominal fit, the x and 
y parameters are fixed to the measured values, and the 
remaining parameter f / 2 is determined to be: 



0/2 



0.00 ±0.02, 



(20) 



where the error is statistical. Figure U shows the angu- 
lar distributions of background-subtracted data and the 
fitted results for \c2 — > 77 events. It is found that all an- 
gular distributions are consistent with the fitted results 
within errors. As mentioned in Section [IV! for the mea- 
surements of the branching fractions we use the formulae 
in Eq. ([5]) to generate MC events for efficiency determi- 
nation of ip' —J- 7Xc2, Xc2 — > 77 decay, with the x, y and 
fo/2 parameters fixed at their measured central values 
\x = 1.55, y = 2.10 and f Q/2 = 0.0). 



In order to test the reliability of the fit, we allow the 
x and y parameters to float in the fit, in which case the 
likelihood fit to data yields 



1.76 ± 0.25, y = 2.00 ± 0.23, 



(21) 



where the errors are statistical. The results are con- 
sistent with the previous BESIII measurements of the 
"0' -» 7Xc2, Xc2 -> 7r+7r~ /K + K~ decays [2(J. 

The goodness of the fit is estimated using Pearson's \ 2 
test [32|]. The data and MC are divided into 6 bins of 
equal width in each dimension (i.e. cos0i, cos 02, 02 ) of 
the three-dimension angular distribution. The numbers 
of events in each cell for data and the normalized MC 
sample are compared. The \ 2 is defined as: 



x 2 = 



E 



(nr_nf C ) 2 

2 



(22) 



where 



,DT 



(cr n DT) is the observed number (its error) of 



signal events after background subtraction in the ith bin 



from data and 



is the expected number of events 



predicted from MC in the ith bin using f / 2 fixed to 
the value determined in the analysis. If the number of 
events in a bin is less than 5, we add the events to the 
adjacent bin. The result of the % 2 test of the fitting is: 
X 2 /n.d.f = 87.3/88 = 0.99, where n.d.f. is the number of 
degrees of freedom. The result shows that the fit quality 
is acceptable. 

Since /0/2 is a ratio, many systematic errors cancel, 
and only the effects due to MC simulation of detector re- 
sponse, the uncertainties on the measured x and y param- 
eters, background substraction, XcO contamination are 
considered here. Among these sources of the systematic 
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-0.6 -0.4 -0.1 ' -0.6 0.2 0.4 ' 0.6 
COS8, 




FIG. 4: Distributions of cos 61, cos #2 and 02 for the 
i>' 7Xc2, Xc2 — ► 77, where the dots with error bars are 
background-subtracted data and the histograms are the fit- 
ted results. 



uncertainties, the MC simulation of detector response is 
dominant; the others are tiny and are neglected. 

As discussed above, the x and y parameters are fixed 
to the measured values from Ref. [26j in the ML fit to 
Xc2 ~ > 77 events in order to obtain the ratio /o/2- In 
the fit we change the x and y central values by one 
standard deviation of the measured values (26j . and find 
that the effect on / /2 is negligible. To estimate the 
uncertainty due to background subtraction, we vary the 
sideband region from (0.07, 0.08) GeV (lower sideband) 
and (0.16,0.20) GeV (higher sideband) to (0.07,0.09) 
GeV and (0.15, 0.20) GeV. After subtraction of the back- 
ground based on the sum of recalculated logarithmic like- 
lihood values, ln£{,, we find that the fitted / / 2 value is 
almost unchanged. From MC simulation, 0.028% of the 
Xco 77 events are distributed under the % c2 signal 
region; the uncertainty due to \c0 contamination is esti- 
mated to be negligible. 

The uncertainty due to the inconsistency between data 
and MC simulation on the angular distributions for \c2 
events can be tested using Xco events. Since the Xco is 
pure helicity-zero, the x and y parameters are expected 
to be zero. In Xco — > 77 decay, the difference of hclicity 
values of the two photons is also expected to be zero, 
so only the helicity-zero term in Eq. ([SJ remains, which 
modifies Eq. ((SJ) to: 



W O (01, 02,02) 



^?/ 2 (l + cos 2 0i ) sin 4 2 + 3a; 2 sin 2 0i sin 2 20 2 - ^^sj/sin26»i sin 2 2 sin202 cos0 2 



\/Hx sin 20i sin 202 cos <j>2 + V&y sin 2 0i sin 2 02 cos 202 + (1 + cos 2 0i) 
"l 



+ /2/0 



V2 



-y' 2 (l + cos 2 0i)(l + 6 cos 2 2 + cos 4 2 ) + 2x 2 sin 2 0i(l + cos 2 2 )sin 2 2 + — icysin20i sin20 2 (3 + cos 2 2 )cos0 2 



~^~ x s i n 201 sin 2 02 sin 202 cos 02 + ~^~V s ' n2 ^1 (•"■ — cos 4 02 ) cos 202 + — (1 + cos 2 0i ) sin 4 02 



(23) 



r 



where the product factor / / 2 is moved to the front fac- 
tor of the helicity-two term and renamed as /2/0, an d 
the (3 cos 2 62 — l) 2 term associated with A = ampli- 
tude in Eq. (J5J) is replaced by 1, so that one can obtain 
the expected angular distribution Wq = 1 + cos 2 d\ from 
Eq. (|23|) if the parameters x = 0, y = and / 2 /o = 0, 
as expected. Therefore, we fit the angular distribution 
°f Xco with the Eq. (|23|) using the same method as in 
Xc2 decays; non-zero x, y and / 2 /o values will indicate 
the inconsistency between data and MC simulation. The 
systematic error is taken as the shift from plus its error. 
The fitted results are x = -0.11 ± 0.09, y = 0.13 ± 0.07 
and /2/0 = 0.00 ± 0.02. The correlation coefficient be- 
tween x and y is -0.27, while it is 0.0 between x (y) and 
fz/Q. Thus we take 0.02 as the systematic error for the 



measurement of / / 2 in the fit to \c2 events. Studies with 
MC simulated data samples demonstrate that a system- 
atic error in modeling the 9%, 2 , and <f>2 efficiency pro- 
duces a shift of approximately the same size for f 2 / in 
Xco sample and / / 2 in \c2 sample, when the latter sam- 
ple is generated with x = 1.55, y = 2.10 and /o/ 2 = 0. 
Therefore, we assume the observed shift from / 2 /o for the 
true Xco data is an estimate of the systematic error on 
the measured values of f / 2 for the two-photon decay of 

Xc2- 
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TABLE IV: The comparison of experimental results for the two-photon partial widths of XcO and \ c 



Quantity 



Si x S 2 x lO^of 

51 x S 2 x 10 5 ( Xc2 ) c 

5 2 x 10 4 (xco) c 
B-2 x 10 4 ( X c2) c 
r 77 ( X co)(keV) 
r 77 (x c2 )(keV) 

n 

fp/2 



PDG global fit results" 



CLEO-c" 



This measurement 



2.16 ±0.18 2.17 ±0.32 ±0.10 2.17 ± 0.17 ± 0.12 

2.24 ±0.17 2.68 ±0.28 ±0.15 2.81 ± 0.17 ± 0.15 

2.23 ±0.17 2.31 ±0.34 ± 0.15 2.24 ± 0.19 ± 0.15 

2.56 ±0.16 3.23 ±0.34 ± 0.24 3.21 ± 0.18 ± 0.22 

2.32 ±0.22 2.36 ±0.35 ±0.22 2.33 ± 0.20 ± 0.22 

0.50 ±0.05 0.66 ±0.07 ±0.06 0.63 ± 0.04 ± 0.06 

0.22 ±0.03 0.28 ±0.05 ±0.04 0.27 ± 0.03 ± 0.03 

0.00 ±0.02 ±0.02 



a The results from the literature have been reevaluated by using the branching fractions and total widths from the PDG 
global fit. 

h The first error is statistical. The second error is systematic error combined in quadrature with the error in the branching 
fractions and widths used. 

c Si = B(ip(2S) -> 7XcO, C 2), S 2 = S(x c o, C 2 -> 77), r 77 (x c o,c2) = r 77 (xco, c2 -> Tf)- 



VI. CONCLUSION 



which the ratio is predicted to be less than 0.5%. 



In summary, we present measurements of the two- 
photon decays of Xc0,2 via the radiative decays ip' — > 
1Xc0.2 ■ We find B( X cO 77) = (2.24±0.19±0.15) x 10~ 4 
and B(xc2 -> 77) = (3.21 ± 0.18 ± 0.22) x 10~ 4 , which 
agree with the results from the CLEO experiment [i~3l |. 
The partial widths r 77 (x c o. C 2) and the ratio 1Z of the 
two-photon partial widths between \ci and XcO are de- 
termined from these measurements. The precision of our 
measurements is improved compared to CLEO's; the fi- 
nal results are listed in Table IIVI 

Since theoretical unknowns cancel in the ratio 1Z, a 
calculation including the first-order radiative corrections 
by Voloshin [H predicts TlW = 0.116 ± 0.010. Our ex- 
perimental result, 1Z — 0.27 ± 0.04, indicates some inad- 
equacy of the first-order radiative corrections that have 
been used to make theoretical predictions for charmo- 
nium decays. 

We also perform a helicity amplitude analysis for the 
decay of ip' — > 7Xc2, Xc2 — > 77; the ratio of the two- 
photon partial widths for the helicity-zero and helicity- 
two components in the decay \c2 11 is determined 
for the first time to be f 0/2 = 0.00 ± 0.02 ± 0.02. The 
helicity-zero component in the \c2 — > 77 decay is highly 
suppressed. This measurement is consistent with the cal- 
culations based on a relativistic potential model [5j, in 
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